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Abstract 

Integral transforms arising from the separable solutions to the Helmholtz differential equation 
are presented. Pairs of these integral transforms are related via Plancherel theorem and, ultimately, 
any of these integral transforms may be calculated using only Fourier transforms. This result is 
used to evaluate the mean value of momenta associated to the symmetries of the reduced wave 
equation. As an explicit example, the orbital angular momenta of plane and elliptic-cylindrical 
waves is presented. 
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I. INTRODUCTION 



The orbital angular momentum (0AM) of light has become a fundamental concept in 



optics 



The fact that light may carry orbital angular momentum which is not re- 



lated to its polarization state has been shown both theoretically [3| and experimentally j4|. 
The optical tweezers toolbox now includes experiments involving the transfer of 0AM from 



paraxial beams js], 0] and vector fields [7[ to dielectric microscopic particles; also, the fea- 
sibility of transfer ring 0AM from laser beams into cold atoms 8|, |9|, and even store it in 
matter ensembles [lO|, has been discussed and shown. Furthermore, a couple of dynamical 
variables of light, involving compositions of linear and angular momenta, have been theo- 
retically identified related to propagation invariant vector fields with transversal elliptic 11 1 
and parabolic symmetries. In the paraxial regime, the elliptic 13| and parabolic 14 1 
beams are well known; they have been generated by diverse means 



a 



15 



angular momentum per photon of the elliptic beams discussed 



17j and the orbital 



18| . Proposed applications 



for these light fields carrying compositions of linear and angular momenta include the optical 



19| and atomic billiards 



20|, 



manipulation of atoms, e.g. thermal atomic cloud splitters 
as well as the optical manipulation of dielectric microscopic particles in elliptic, hyperbolic 
and parabolic trajectories in optical tweezers. The transition from mere optical manipu- 
lation towards optical control of matter, and further uses in other fields, e.g. in quantum 
information protocols 2l|, |22j , implies measuring the linear and orbital angular momenta, or 



compositions of these, of arbitrary light fields. Here, a simple schema to calculate momenta 

brm, based on the Hilbert space model for 



23 



is presented for the benefit of ad- 



for scalar wave fields through the Fourier trans 
separation of variables of Helmholtz equation 
vanced undergraduate and graduate students interested in the 0AM of light and related 
topics. 

In the following section, a brief overview of the Helmholtz equation, its symmetries 
and separable solutions, is given; the separable solutions for Cartesian, circular-cylindrical 
and elliptic-cylindrical coordinate systems are shown as an example. In the third sec- 
tion, wave decompositions and integral transforms related to the Helmholtz equation are 
treated; Plancherel theorem, the keystone of the proposed schema, is presented in this sec- 
tion. Fourier and Circular-cylindrical wave transforms of the separable solutions presented 
in section 2 are calculated for the sake of clarity. Next, the actual schema to calculate 
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momenta through Fourier transforms is discussed and the mean value of the orbital angular 
momentum for Plane and Elliptical waves is calculated. Finally, a summary and plausible 
applications are discussed. 



II. SYMMETRIES OF THE HEMLHOLTZ EQUATION 

The Helmholtz differential equation, 

(v2 + k^)^^J = 0, (1) 

where the symbol is the Laplacian operator and k is the magnitude of the wave vector 
k = fcfsin-j^ cos g9, sin ■i^ sin 09, cos'j^), and ^ and (f are the Euler angles, is separable in eleven 



coordinate systems 25142 7l|. Each of these coordinate system defines a three element set of 



operators that span a Lie algebra 



23 



24| . These operators fulfil the eigenvalue equation 



Ofvl/ = of vl>, (2) 



with z = 1, . . . , 11, and j = 1, 2, 3. By defining a characteristic label set given by the three 
eigenvalues corresponding to a given coordinate system, Aj = it is possible to label 

the set of functions that solve the reduced wave equation in the z-th coordinate system as 
{\E'Ai}- A common element to all characteristic label sets is the square linear momentum 
that fulfils p^^E^Ai = — V^\E'a. = k^^^A^; ^-c, the first element of each and every characteristic 
label set A^ is the magnitude of the wave vector, of = k. The four coordinate systems 
with axial symmetry, Cartesian (i = 1), circular- (i = 2), elliptic- (i = 3), parabolic- {i = 4) 
cylindrical, share the 2;-component of the linear momentum as the second dynamical variable, 
Pz'^Ai = {k COS 6)"^ j>^.; i.e., for these four coordinate systems the second element of each of 
the characteristic label set is the Euler angle 6, of = 6 with i = 1, ... ,4. There are five 
coordinate systems that share the orbital angular momentum, Iz'^Ai = n'^Ai, with 1 = r x p; 
for these, the second element of each of the characteristic label set is the discrete orbital 
angular momentum number n, also known as topological charge. Finally, there are another 



thirteen operators that can be easily found in the literature [23|, |24 |. 

In optics, there is a particular interest in propagation invariant scalar wave fields. There 
exist four fundamental separable families with this characteristic; for the sake of brevity, we 
will present just the cases for plane, circular- and elliptic-cylindrical waves. 



3 



A. Plane waves 



In the Cartesian coordinate system, x = (x,y,z), the separable solution of the Helmholtz 
differential equation are plane waves, 

^Ai = (Sin?9)^/^ grfc(xsin^cosv.+!/sintfsinv+^cos^)^ 0<^ <n, 

(3) 

The conserved operators are the linear momenta Py and Pz, 

Py^A, = (A;sin^?sin(/?)^fAi, (4) 

P.^'ai = (A;cosi9)«'Ai, (5) 

yielding a characteristic label set given by the magnitude of the wave vector and the two 
Euler angles, Ai = {k,^,(p}. 

B. Circular- cylindrical waves 

For circular-cylindrical coordinates, x = (r cos0, r sin 0, z), the separable solution of the 
Helmholtz differential equation are circular-cylinder waves, 

^A2 = ^''(27rsin^9)^/^ J„(A:sin^9 r)e*("'^+'=""°''^\ < ^9 < vr, n = 0, ±1, ±2, . . . 

(6) 

where the function Jn{-) is the Bessel function of the first kind. These scalar wave fields, so 
called Bessel waves, conserve the linear momenta p^ and the orbital angular momenta 1^, 

P.^Aa = (A;cos?9)^A2, (7) 
k^A, = n^'Aa, (8) 

thus the characteristic label set is given by the magnitude of the wave vector, one of the Euler 
angles and the orbital angular momentum number, or topological charge, A2 = {k,^,n}, in 
that order. 

C. Elliptic-cylindrical waves 

In elliptic-cylindrical coordinates (x = / cosh ^ cos 77, y = /sinh^sinr/, z), the separable 
solution of the Helmholtz wave equation are the so called even and odd Mathieu waves, in 
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that order, 

^ca3 = (sin^9)^/^c„Ce(e,g)ce(r/,g)e*'=^™^'', < ^9 < tt, n = 0, 2, . . . (9) 
^SA3 = (sin^?)^/^s„Se(e,g)se(r/,g)e*'^''=°''', < < vr, n = 1, 2, . . . (10) 

where the normahzation coefficients are given by the expressions 

•^2" = T^'^^2n(0;g)ce2n(vr/2;g), (11) 



^0 



C2n+i = 45;^ce2„+i(0;g) 9„ce2n+i(t^;g)L=^/2, (12) 

S2n+i = j^;^ (9„se2„+i(t;;g)|„^oSe2„+i(7r/2;g), (13) 

^2n+2 = — 7^;^ 9„Se2n+2(t';g)L=0 ^^Se2n+2(f;g)L=^/2, (14) 
?^2 

where ce„, (se„) and Ce„, (S'e„) are the ordinary and extraordinary n-th order Mathieu 
functions of even (odd) parity, sometimes named angular and radial Mathieu functions, 
in that order. The separation constant is given by g = (//c sin^9/2)^ with 2/ being the 
inter-foci distance that defines a particular elliptic cylindrical coordinate system. Mathieu 
waves conserve, apart from the 2;-component of the linear momenta, a linear superposition 
of orbital angular momenta Iz and linear angular momenta px 

[ll + /^p^.)^CA3 = a^vl/cf,),^ (15) 
{ll + Pp>s^, = (16) 

such that the characteristic label set is given by the magnitude of the wave vector, one Euler 
angles related to the 2;— component of the wave vector and the discrete order number, A3 = 
{k^d^n] in that order. The even (odd) eigenvalue a„ (6„) is fixed for a given combination 
of n and q. 



III. WAVE DECOMPOSITIONS AND INTEGRAL TRANSFORMS 

Each one of the eleven sets of separable functions, {\1/a,}, that solve the reduced wave 
equation forms an orthonormal basis, 

j rfx ^r^^(x)vl/^,(x) ^ (M^aJ^a;) = m - AO, (17) 



where the bracket, (■, ■), notation is used as a shorthand for integration over configuration 
space, the upper bar notation, ^, stands for complex conjugation and the notation 6{-) stands 
for a combination of Dirac and Kronecker dehas particular to each coordinate system. A 
'$j>^.-w8ive decomposition of an arbitrary well behaved function $(x) can be defined for each 
of the eleven orthonormal basis, 



$(x) = J rfA, 0(A,)*A,(x), (18) 

alongside an integral transform, 

0(A,) = y"rfx ^^A,(x)<|.(x) = (19) 

It is possible to show that each and every one of these eleven wave decompositions can 
be related to each other through their amplitudes, by substituting the corresponding wave 
decomposition, Eq. (|T8l) . of ^^Ai(x) and $(x) into the integral transform given by Eq.([T9]). 



0(A,) = J rfA,^A,(A,)0(A,.), (20) 

In shorthand, 

(v1/aJ$) = ((V^aJ0))a,, (21) 

where double angle-brackets, {{■\-))Aj, have been used to denote integration over the phase 
space defined by the characteristic label set Aj. The above expression, Eq. (!2T]) . is the well 



known Plancherel theorem 



28 



30l|: it is particularly useful when both the integral transforms 
related to a \l/Aj-wave decomposition of $ and \I^Ai are known. 

Unitarity of the eleven Helmholtz-related integral transforms follows straightforward from 
Eq. fl2Tl) yielding Parceval theorem, 

($1$) = ((0|0))a.. (22) 

A. Plane wave decomposition and Fourier transform 

Take as an example the plane wave decomposition 

$(x) = J c/Ai 0(Ai)(sin^?)i/2g*k.x^ (23) 
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related to the Fourier transform 

4>(k) = J (sin^9)^/2g-rk.x^(^^)_ ^24) 

The Fourier transform of a plane wave is given by the expression 

V^a;(Ai) = {sm^)-'^'S{^-^')S{'d-'d'), (25) 

where the function 6{-) is the Dirac delta function. A circular-cylindrical wave has a Fourier 
transform 

V'Ai,(Ai) = (27rsin^9')"'^'e*"'^<5(e-79')- (26) 
For elliptic-cylindrical waves, Fourier transforms are 

V'CA.(Ai) = {7ism^)-'/'cen{^,q)6{^-^'), (27) 
V'SA^(Ai) = {7rsm^)-'^^sen{^,q)6{^-^'). (28) 

B. Circular-cylindrical wave decomposition and related integral transform 

Now, take the circular cylindrical wave decomposition, also known as Bessel wave decom- 
position, 

/•Tr 

$(x) = V / sin^?rf^9 0(A2)«"(27rsini9)i/V„(A;sini9 p)e*("^+'="^°'''). (29) 

n 

This circular wave decomposition defines a Bessel transform related to the Helmholtz equa- 
tion; let us call it a Helmholtz-Bessel transform (not to be confused with the Fourier-Bessel, 
so called Hankel, transform) for the sake of having a short name for it, 

0(A2) = (^aJ$) = y"rfx $(x)(-z)"(27rsin?9)^/V„(A;sin^? p)e-*("^+'="^°'^''). 

(30) 

Sometimes, it can be bothersome to calculate the value of this expression. If the Fourier 
transform </'(Ai) of the function $(x) is known, it may be simpler to calculate the coefficients 
of the Bessel wave decomposition, by using Plancherel theorem, Eq.f l2ip . as 

(/)(A2) = (27r)"^/2^sin?9)^/2 / (t){k,d,^)e-''''<' . (31) 



For example, the Helmholtz-Bessel transform of a plane wave, described by the label set A'^^, 
is given by the expression, 

7/^A',(A2) = {271 sm^y^/\-'''^5{'& -'&'), (32) 

while that of a Bessel wave is given by, 

V'A,(A2) = (sini9)-^/'5(^-i9')5n,n', (33) 

where the symbol 6.^. stands for Kronecker delta. And that of a Mathieu wave is 

^CA.(A2) = 2"l/M("J^(„,^„,2)(g')5(^-^')^n,2^+Kn.od2), (34) 
^.A.(A2) = 2-l/2i?(':;)_(„,^„,2)(?')5(^-^')^n,2,n+Kmod2), (35) 

with m = 0, 1, 2, . . . and the coefficients A^*^ ( -B]*'') are the j-th coefficients of the cosine 
(sine) expansion of the ordinary Mathieu functions of i-th order. 

IV. QUANTIFYING MOMENTA THROUGH THE FOURIER TRANSFORM 

It is possible to calculate the mean value of one of the aforementioned sixteen conserved 
operators Oj''\ as well as linear superpositions of them, for any given wave \1/ using its 
\I/Ai-Wcive decomposition, 

1$) = J]o«|(vl/Aj$)r. (36) 

Via Plancherel theorem, Eq. (!2Tl) . it is possible to give this expression in terms of Fourier 
transforms, which may be already available and provide an easier path to calculate the 
overlaps, 

1$) = 5^o«|((^aJ0))a,P, (37) 

A, 

where the phase space defined by the label set Ai = {k, {}, (p} yields 

roc nix 

m<P))M= dk sini?di?^A,(Ai)0(Ai). (38) 

Jo J-TT Jo 
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A. Orbital angular momentum of plane waves 

It is simple to calculate the well-known null mean value of the orbital angular momentum 
for a plane wave, Ai = {k, ^, (p}, by using Eq. fl57|) and Eq. flHT]) . 

oo 

(v1/aJ/,|x[/^,) = 5^ n|(fel^Ai))A'J', (39) 



?i=— oo 
oo 



= Yl ^|(sin^9)'/'(27rsini9)-^/2g«n,?^^^_^^|2 ^^q^ 

n=— oo 

= 0. (41) 
B. Orbital angular momentum of elliptic-cylindrical waves 

Again, by use of Eq. fl37|) and Eq.f l34|) . it is possible to calculate the mean value of the 
angular momenta for any given elliptic-cylindrical wave as 

oo 

(vl/cA3|/,|vI/cA3) = 5^ n |((V;aJ^ca3))a,|2, (42) 

n=— oo 

°° 1/2 

= J][m+(nmod2)/2]| (sin^/sin^) 4':!,+(nmod2)^(^ - 

oo 

= ^ [m + (nmod2)/2] l4:l+(n.od2)^(0)P, (43) 

m=0 

oo 

(vl/.A3|/,|M/SA3) = l((^Aj^SA3))A;r, (44) 



n=— oo 
oo 



= ^ [m+ (nmod2)/2] | (sin ^/ sin ^) i?tV{nmod2)^(^ " 

m=0 
oo 

= 5^ [m + (nmod2)/2] |5£+(_,,)5(0)r. (45) 

m=0 

Note that the divergences in all these expressions are related to the fact that ideal scalar wave 
solutions to the differential Helmholtz equation have a divergent normalization, Eq.( ll7p . 

V. FINAL REMARKS 

A schema to calculate momenta related to the symmetries of the Helmholtz differential 
equation has been presented. The procedure relies on the orthonormal bases and the sixteen 
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conserved operators derived from the eleven coordinate systems where the reduced wave 
equation accepts a separable solution. It has been shown that all the integral transforms, 
or wave overlaps, yielding the coefficients for the wave decompositions related to the afore- 
mentioned orthonormal basis can be calculated through Fourier transforms via Plancherel 
theorem. These overlaps are particularly useful to calculate the mean values of the con- 
served operators, and linear superposition of them, constructed from the three linear and 
three angular momenta. As an example, the mean value of the orbital angular momentum 
of Elliptic-cylindrical waves, so called Mathieu waves, has been presented. With further 
efforts the schema can be extended from these ideal scalar wave fields to real world beams; 
jy using the orthonormal bases given by Gauss, Bessel-, Mathieu- 31| and Weber-Gauss 



32| beams it may be possible to implement computational procedures to calculate momenta 



of captured experimental images. An extension of these procedures to non-paraxial optical 
vector fields is plausible. Moreover, integral transforms and decompositions are of use for 
those in the quantum optics fields as they define decompositions between boson modes; e.g. 
expressing creation/annihilation operators of Bessel modes in terms of creation/annihilation 
operators of Plane wave modes 33|. 

Furthermore, to give just an example from outside the field of optics, in quantum me- 
chanics, the treatment of molecular torsional levels for water and similar molecules is related 
to Ince polynomials that are a solution to the Helmholtz differential equation in paraboloidal 



coordinates 34] • In these cases, the presented schema can help in calculating suitable com- 



binations of linear and orbital angular momenta of the different torsional eigenstates of the 
molecules. 
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